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Abstract 

The two dominant components of the cosmic budget today, pressureles matter and dark energy, 
may or may not be interacting with each other. Currently, both possibihties appear compatible 
with observational data. We propose several criteria based on the history of the Hubble factor that 
can help discern whether they are interacting and whether dark energy is phantom or quintessence 
in nature. 
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I. INTRODUCTION 



Recently, Sahni et al. used the luminosity distance di = c(l + z) dz'/H{z'), valid in 
spatially fiat Friedmann-Robertson- Walker (FRW) universes, to propose a diagnostic to tell 
the cosmological constant from evolving dark energy fields For universes dominated by 
cold matter (CM), subscript m, and dark energy (DE), subscript d, with constant equation 
of state, w = Pd/ Pd < —1/3, Friedmann's equation takes the form 

E\x) = nmox' + ndox"^''-'"\ (1) 

where x := 1 + z = ao/a, and E{x) := H{x)/Hq is the normalized Hubble function. In 
writing this expression it was implicitly assumed that cold dark matter (CDM) and DE 
evolved separately. As usual, f2jo := pjo/(3-f^o) with i = m,d stand for the energy 
densities of the components in units of the critical density, := 8ttG , and a zero subscript 
denotes the present-day value of the corresponding quantity. 
Sahni et al. [ij introduced the redshift dependent function 



Omix) := ^^^^ , (2) 



which for spatially fiat FRW cosmologies {Qmo + ^do = 1) reduces to Om{x) = Qmo + 
(1 — ^lmo){x^^'^^'^^ — l){x^ — 1)^^. It has the interesting feature that for the ACDM model 
{pd = constant, w = —1) , it yields Om(x) = Qmo , whereas for quintessence fields {—l<w< 
— 1/3) and phantom fields {w < —1) it gives Om{x) > flmo and Om{x) < Qmo, respectively. 

Since this diagnostic does not depend on the present values of the density parameters, fi^o 
and Qdo, it circumvents the drawbacks brought about by our comparatively poor knowledge 
of them. In fact, it is expected to be useful when accurate data of H at different redshifts 
become available. This seems more feasible than to substantially upgrade our knowledge of 
flmo and fldo in the near future. 

At present, the ACDM model fits the observational data, within statistical errors, rather 
well -see, e.g. {2]. Nonetheless if future data come to reveal that Om{x) ^ Qmo , a number 
of DE alternatives should be considered including scenarios in which the DE interacts non- 
gravitationally with other components. 

Here, we assume that DE and CDM, subscript c, interact with each other (whence they 
do not evolve separately) but not with baryonic matter, subscript b, since interactions with 
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y con- 
61 and 



ba^yons are strongly constrained by experimerrts g This possibility ,s bemg aetive 

sidered in the literature -see and references therein- since it is rather natural [5|, 
more general than otherwise, it substantially alleviates the coincidence problem (i.e., "why 
are the densities of matter and dark energy of the same order precisely today"), and agrees 
with observation. 

In this paper we propose several criteria, based on the evolution of the Hubble factor, to 
help discriminate models in which DE and CDM interact with each other non-gravitationally 
from models in which they do not, and to tell quintessence DE energy models from phantom 
models. Clearly, for these criteria to be useful accurate data of H{z) are needed. While 
the present data exhibit large error bars 7H9| the situation may well improve comparatively 



soon thanks to upcoming and planned experiments. The Baryon Oscillation Spectroscopy 
Survey (BOSS) 5-years project 1^ aims at measuring the absolute cosmic distance scale 
and expansion rate with percent-level precision at redshifts z < 0.7 and z ~ 2.5 using 
the standard rule furnished by the baryon acoustic oscillations. The forecast precision for 
H{z) a.t z = 0.35, 0.6 and 2.5 is 1.8%, 1.7% and 1.2%, respectively. Data at intermediate 
and larger redshifts will partly come by adapting to 2; > the method proposed by Hogan 



ll| which makes use of the gravitational radiation emitted by black hole binary inspiral 



events to determine Hq with a precision better than 1%. The method may be extended 
up to z ~ 1 redshifts, and even to 2; ~ 10 by observing massive black hole binaries since 
there will likely be a direct electromagnetic count^erpart of the merger [l^. An even more 
promising method, proposed by Corasaniti et al. [13] to constrain cosmological parameters, 
exploits the Sandage-Loeb test [14]. It consists in measuring the redshift of distant objects 
at two separate times, say to and to + Atg, at the present epoch. The redshift variation 
of the source (subscript s) between these two instants is, at a very good precision, A^;^ ~ 
(d(to) — ci{ts))/ci{ts) -see Eq. (5) in [l3^. This expression, re-written as H{zs) — {Azg/ Ato) — 
Ho (1 + Zs), provides us with the value of the Hubble factor at the time the source emitted 
the electromagnetic signal that is currently arriving at our detectors (either telescopes or 



radio-telescopes), modulo we know Ho accurately enough. The authors of [13[ suggested 
using the absorption lines of the Lyman-a forest in the redshift interval 2 < z < 5 and 
Ato = 10 years. 

Altogether, encouraged by the prospects of a greatly improved knowledge of the history of 
H{z) we present several criteria that may help telling apart interacting from noninteracting 
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dark energy models. 

Section II presents on phenomenological basis some possible expressions of the interaction 
term. Sections III and IV introduce various criteria to distinguish interacting from nonin- 
teracting DE models. Section V focus on the statefinder parameters and compares their 
usefulness with the criteria based on the Hubble history. Finally, section VI summarizes our 
findings. 



II. EXPRESSIONS FOR THE INTERACTION BETWEEN DARK ENERGY AND 
DARK MATTER 

Since the dark components are assumed to interact with one another also non- 
gravitationally but not with baryons the conservation equations take the form 

Pb + 3i/pfe = 0, (3) 
Pc + 3Hp, = Q, (4) 
Pd + 3Hil + w)pd = -Q, (5) 

where Q denotes the interaction term and the relationship pb + Pc = pm is understood. It is 
worth noticing that when Q > -as we shall consider throughout- the coincidence problem 
is alleviated since pc decreases more slowly with expansion and pd more quickly, and the 
problem may even be solved in full jl^. On the other hand, if Q were negative, the second 
law of thermodynamics could be violated [l6], and the energy densities of CDM or DE could 
become negative at high or low redshifts. 

In the absence of a fundamental theory for dark energy the quantity Q cannot be derived 
from first principles. However, we may guess likely expressions for it by noting that Q must 
be small [17| (at least lower than 3Hpm) and depend on the energy densities multiplied by 
a quantity with units of inverse of time, a rate. We shall deal with two different possibilities 
considered in the literature, namely, (i) that the said rate is proportional to the Hubble 
factor, and (ii) that it is just a constant jisl . 

In case (i), we have that Q = Q{Hpc, Hpd). By power law expanding this function and 
retaining just the first term one follows Q — Cc Hpc + Hpd- Given the lack of information 
about the coupling, it appears advisable to work with just one parameter rather than two. 
Thus, three different choices arise, namely, ec = 0, = 0, and ec = e^. It is worth 



mentioning that expressions of this type can be obtained from the scalar-tensor theory of 



gravity developed by Kaloper et al. [l9| -see e.g. 20|, |2l|. In the next section we shall 
consider successively three expressions for Q with < e << 1, while keeping w = constant 
throughout, and look for diagnostics -based in the history of H{z) that may tell interacting 
from noninteracting models. 



In case (ii) we have Q = Q{pc, Pd) and proceeding as before we write 

(5 = 3 {TcPc + ^dpd) , 



(6) 



c, d) quantities denote constant rates. This expression was inspired in 



where the Fj {i - 

reheating models 22|. curvaton decay 23|, and decay of DM into radiation 2J]. Regrettably, 
one of the energy densities (pm or pa) becomes negative either at early or late times [l8| and, 
generally, it does not solve the coincidence problem. Nevertheless, we will briefly consider 
it in the section after next. 



III. INTERACTION TERM Q = QiHpc, Hpa) 

We begin by assuming that the interaction term takes the form 

Q = SeHpc , 

where the factor 3 was introduced for mathematical convenience. 



(7) 



A. Decaying cosmological constant 

One way to explain the present small value of the cosmological constant is to assume a 
transfer of energy from the quantum vacuum to particles and/or radiation -i.e., a decaying 
cosmological constant (DCC). Along the years a diversit y o f phenomenological models have 
been proposed, many of which are summarized in Ref. 25|]. If the quantum vacuum {w = 
— 1) slowly decays into non-baryonic CDM particles, integration of Eqs. ([3])- ([5]) yields 



pb 


= Pbox^, 


pc 


= PcQX^^^ 


Pa 


= PA0+ \ 



PcO [X 



3(l-e) 



1) 



(8) 
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This kind of models are not purely phenomeno 
some persuasive backings from the theoretical side 



ogical any more for they have received 



26 



-30l. which makes their consideration 



all the more interesting. In particular, by means of the renormalization group theory the 



authors of 27|] arrive to a simple expression for e {u in their notation) in terms of the masses 
of the super- heavy sub-planckian particles -see their Eq. (4.10). A rather similar expression 
was found by removing the flat space-time contribution to the vacuum energy of quantum 
fields in homogeneous and isotropic universes by employing a technique similar to that used 



in the computation of the Arnowitt-Deser-Misner mass -see Eq. (39) in Ref. j29|. Moreover, 



consistency both with primordial nucleosynthesis and large scale structure formation leads 
to the constraint | e |<^ 1 jsil- Recently, DCC models were revisited by Basilakos 



et al. 32| and, for some specific choices of A(t), confronted with observational data from 
supernovae type la, baryon acoustic oscillations, the cosmic microwave background shift 
parameter, and the growth rate of galaxy clustering. 

As said above, the coincidence problem gets alleviated since 



|r/r 



DCC 



\r/r\ 



1 - e + ro] < 1 



\ACDM / 

i.e., currently the ratio r = {pf, + Pc)/pa varies more slowly than in the ACDM model. 
This also holds true for the other interaction terms considered below. 
Using Friedmann's equation, 3if^ = [pb + Pc + Pa] , we can write 

E\x) = Qhox' + Qcox^^'-'^ + ^Ao + - l] , (9) 

where ^ao = 1 — ^bo — ^co- Obviously, in this case Om{x) ^ fimo with f2mo = + ^co, 
but rather 



Om{x) = n^o + ^ ^3 , . (10) 



1 — e — 1 

Thus, the Om(a;) criterion would mistake a DCC model by a quintessence field (phantom 
field) if applied at low redshifts (high redshifts). Similar situations occur if the said criterion 
is used on models considered below in this paper. 

The first and second derivatives of E'^{x) with respect to are 
dE'^{x)/dx^ = ilbo + ^cox'^", and <PE'^{x) / d{x^Y = —eVlcQX~^^^~^''\ respectively. Therefore, 
for large redshifts the first derivative tends to Vtbo from above. Likewise, E'^{x) is concave 
(i.e., d"^ E"^ (x) / d{x^Y < 0). By contrast, for noninteracting dark energy fields with w = 
constant -see Eq. ([1])- one has, dE'^{x)/dx^ = flyno + (1 ~ ^mo) (1 + 'w)x^^. For large 



redshifts this expression tends to Qmo from above if DE is of quintessence type and from 
below if it is phantom. In view that, observationally, Qbo and Qmo are separated by a non 
small gap (the latter might be about six or seven times larger than the former) this may 
serve discriminate DCC models from noninteracting quintessence and phantom models. One 
may think that we would have arrived to the same result just just by dividing E'^{x) by x^; 
i.e., E^{x)/x^ = Qmo + (1 — ^mo) x^^. However, this is not true because of the factor (l + w) 
in the second term of the derivative which makes the latter tend to Qmo faster, and it is key 
to tell whether the DE is quintessence's or phantom's type. 

As a corollary, we may say that the sign of the second derivative, may discriminate DCC 
from noninteracting phantom models, the former being negative and the latter positive, 
d'^ E"^ (x) / d{x^y = (1 — fi^o) ""^(l + u)) x^^^~^^] only that this criterion can be used also at 
smaller redshift than the preceding one. 

Likewise, by measuring H{z) at two separated redshifts, zi and Z2 with zi < Z2, they 
may be discriminated. Indeed, for noninteracting phantom models the quantity A(xi, X2) := 
E^{x2) - E^{xi) is 

A(a;i,X2) = fi™o [xU^ - x^) - )] + - (11) 

and it can be negative if the redshifts are sufficiently apart. By contrast, the corresponding 
quantity for DCC models, 

A(xi, X2) = Qbo {xl - xl) + Qco (y^) [^2^'"'^ - ^1^'"'^] ' (12) 

will always be positive-definite. 

Let us assume that, thanks to extra information, we come to know that the right model 
must be a DCC one. Note that to discriminate which one (i.e., its e value) neither direct 
use of last expression nor of Eq. will prove very useful, given the comparatively wide 
uncertainty in the value of Qco- However, the following strategy may work. Consider two 
competing DCC models (say A and B) and recall that Qf,Q is reasonably well determined 
from primordial nucleosynthesis (0.036 < Qbo < 0.047, see, e.g. Ref. [33]) together with the 
present value of the Hubble factor. The ratio 

_ [A{x,,X2) - n^oixl-xl)]^ -f ^-'A 4'"'^ - xf-'^^ . . 

^ ■ [A(xi,X2) - {xl - xl)]^ \l-eAj xf"^^ - xf^^^ ^ ^ 

depends on ca and €3 only. This will enable us to determine which DCC model, A or B, fits 
better the H{z) data. 



It should be noted that, in this regard, it looks more advantageous to use this other ratio 

^ Om{x2) - ^bo ^ 3^1-1 xf^''^ - 1 
""^"^ - Om(xO -nbo~ xl-1 - 1 ' ^ ^ 

since it can constrain e more easily. 

At any rate, aside from this specific instance, the Om{x) function (Eq. ([2])) is not much 
useful when dealing with interacting models since, in these models, CDM and baryonic 
matter redshift at different rates with expansion whereby the current fractional densities, 
QbQ and Qco, enter the expression for E"^ multiplied by different powers of x. This explains 
the restricted application range of the said function and we will not resort to it any more in 
this paper. 



1 . Decaying dark energy field 

For the more general case w = constant 7^—1, one has 

Pb = Pbo x^ , 

Pc = Pco x^^^~^^ , 

P, = P,ox'^'^-^ + Pco [x'^'^-^ - . (15) 



Thereby 



E'{x) = n^ox' + n^ox'^'-'^ + n,ox'^'^-^ + (^-] n^o [x'^'^'^^ - x'^'-^^] , (i6) 

\e + w / 



being 



' ' Qbo + ^co{l - e)x-''' + Qdoil + w)x^^ 



dx^ 



+ (^-) [{I + w)x'^ - (1 - e)x-^^-] , (17) 



and 



£E' 



d{x^Y 

+ ( ^— ^ ^],o k(l + w)x''^'"-^^ + e(l - e)x-=^(^+^n . (18) 
\e + w J 

In view of the hypothesis < e << 1 and the fact that observation reveals that << 
1, the right hand side of (fT7j) is dominated by its two first terms. Likewise, the first term 
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dominates the right hand side of (fT8l) . Therefore, for not large redshifts (i.e., z < 3) the 
normahzed Hubble function, E'^{x), is growing and concave. In consequence, its shape will 
not tell a decaying energy field (with w ^ —1) from DCC models or from noninteracting 
quintessence or phantom models with constant equation of state. However, discrimination 
will likely come at higher redshifts since, as before, the first derivative may tell interacting 
from noninteracting cosmologies depending on whether dE'^{x)/dx^ tends to Qbo or Qco at 
large redshifts. 

Figure [1] shows the evolution of the second derivative of E'^{x), as given by Eq. (fTSl) . for 



the best fit values of w and e of model HI of He et al. [3J] (dotted line) as constrained by 
data from supernovae, baryon acoustic oscillations, cosmic microwave background, and the 
value of the Hubble constant, Hq. Also shown is the noninteracting case (solid line). The 
fact that both lines nearly overlap emphasizes the need for accurate data of H{z). 



0.00 
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-0.15- 



-0.20 - 
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FIG. 1: Dotted line: plot of the second derivative of E"^ in terms of for the best fit values, 
w = -1.02, e = 0.0006, with Q^o = 0.230, and n^o = 0.725 of model III, Q = 3Hepc, of Ref. 
Solid line: same as before but with e = (no interaction). 
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In the present case 



X 



3(l-e) 



— X 



3(l-e) 
1 



cO 



X 



3(1+M«) 



— X 



3(1+10) 
1 



— X 



X. 



3(l+ii)) 



— X 



3(1+10) 
1 



I 3(l-£) 



(19) 



e + w ^ 

For < e << 1 and 1 + w ~ the two last terms are sub dominants, and the ratio defined 
by the first equahty in ( 1T3|) comes to be 



3(1-£a) 3(1-£a) 



X 



3(l-€fl) 

■^2 



3(l-eB) ■ 
Xi 



(20) 



Thus, thanks to the factor (1 — eB)/(l — ca) on the right hand side of Eq. f fT3l) we hope 
discriminate DCC models from interacting dark energy models with w close to the cosmo- 
logical constant value, —1, and also discriminate within models featuring the interaction (|7]). 



B. Interaction term proportional to H 

We now assume that the interaction term takes the form 

Q = SeHpd (21) 



which may be key to explain the non- vanishing temperature of sterile neutrinos 35|, 36 |. 
Then, the conservation equations (EJ-dlj) integrate to 



Pb = Pbox^ , 



Pc = Pcox=^+ f^— ) Pdo x\ 

Pd = Pdox=^(i+'"+^). (22) 



Hence, 



E^(x) = x^ -n 



(^) - • (23) 

For w = —1 (i.e., DCC with interaction given by Eq. fl2Tl) ). it reduces to 

E\x)=x' - ^^[x^-x^^] , (24) 

which is to be compared with (Q. Again, E'^{x) is growing and concave whence its graph 
will tell neither DCC models with interaction term given by ( 12T|) from DCC models with 

10 



interaction term given by ([7]) nor these models from noninteracting quintessence models with 
w = constant. 

Expression fl24p will help constrain e without need of knowing ^ao since the ratio 

[iz/ yx) X \x=xi _ Xi Xi , , 

\F'^('r\ — t31 ~ _ „3e \-^^) 

[J^ l^xj A ix=X2 •'^2 "^2 

does not depend on that quantity. 

Unlike the two previous cases, as inspection of the right hand side of Eq. ( I23l) shows, 
the behavior of the derivative dE'^{x) / dx^ at high redshifts will not discriminate interacting 
from noninteracting cosmological models. For the discrimination to be possible e should 
take unrealistic high values. 

By deriving twice ( 123|1 we get 

'^'^'^''^ ndow{l + w + e)x^^^^'-^\ (26) 



d{x^f 

When this latter quantity gets accurately determined, we will be able to set useful con- 
straints on w and e since its ratio for two decaying dark energy (DDE) models (say A and 
B) does not depend on the VL parameters, 

d'^E^dix^fU ^ WA{l + WA + tA) 

Once again E^{x), given by ( l23l) . is a growing function of redshift, concave for quintessence 
fields and convex for phantom fields. Accordingly, the shape of E'^{x) can tell decaying 
quintessence fields from decaying phantom fields but not from noninteracting DE fields. 

Figure [2] depicts the evolution of the second derivative of E'^{x), as given by Eq. (l26l) . 



for the best fit values of w and e of model II of Ref. [3J| (dotted line) as constrained by 
observational data. The solid line corresponds to the noninteracting case (solid line). As it 
is apparent the lower the redshift, the easier the discrimination. 

C. Interaction term proportional to H {pc + Pd) 

Assuming the interaction term 

Q = SeHip, + Pd) , (28) 
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0.016- 

d'E'/d(x'f 



0.008 - 




FIG. 2: Dotted line: plot of the second derivative of E'^ in terms of for the best fit values, 

—I 

w = -1.03, e = 0.024, with Qao = 0.717 of model II, Q = SHepd, of Ref. [M]- Solid line: same as 
before but with e = (no interaction). 



widely considered in the literature -see, e.g. 
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, uM, |37[-, equations ([3])- ([5]) integrate to 



Pb = Pbox^ , 



.72 



respectively. 
Here 



Pd = ^ [-CM + B) + C2{B- A) , (29) 



2 2 

7i = - [2 + w + ^/w{w + 4e)] , 72 = - [2 + ti; - ^/w{w + 4e)] , 



A = {w + 2e), and fi=Vw(w + 4e). 

As in the two previous cases, the three energy densities remain semi-positive definite for all 
X. Therefore, 

E''{x) = Qw +^ [Qc, [2e -{A + B)] - Qc^ - B) x'^^} + Qc^ , (30) 

12 



where 

nc, = ^[nco{B-A)-2endo], and Qc^ = [^co {A + B) + 2 eQao] , 

Z JD Z ID 

also stays non- negative as it should. Further, E"^ is a growing function of x^, concave for 
quintessence and convex for phantom -see Fig. [3]- which will help discriminate phantom from 
quintessence DE in interacting models with Q given by Eq. (l28ll . However, as is apparent, 
for > 5 the graphs show a nearly straight line behavior whereby one must focus on 
redshifts below, say, 1.7. 

Inspection of the right hand side of ( 130|) shows that for interacting models (e > 0) the 
derivative dE'^{x)/dx^ tends to Qbo at high redshifts. Accordingly if observation reveals that 
behavior, it would be suggestive of interaction. 




FIG. 3: Left panel: Graphs of E"^ vs. as given by Eq. (j30p for three values of w. Right panel: 
Graphs of the second derivative of E'^ with respect to for the same three values of w. Notice that 
only phantom models present convex curves for reasonable values of the parameters. In drawing 
the graphs we took = 0.04, QcO = 0.26, ^do = 0.70, and e = 0.05. 



Figure m shows the evolution of E"^ for the best fit value of model IV (e = 0.0006) of [34 1 
(dotted line) and the values corresponding to the la error, e = 0.0011 (dashed) and e ~ 
(solid); obviously the latter practically coincides with the noninteracting case. Likewise, Fig. 
[5] shows the evolution of (P E"^ [x] / d{x^Y , for the best fit values of the same model (dotted 
line) and for the noninteracting case (solid line). As is apparent the graphs practically 
overlap. 
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5 3 6 

X 

FIG. 4: Dotted line: plot of E'^ in terms of for the best fit values, w = -1.03, fico = 0.230, 
and Q.dQ = 0.725, e = 0.0006, (dotted line) and the e values corresponding ^to the la errors, 
0.0011 (dashed) and ~ (soUd) of model IV, Q = 3He{pc + pd), of Ref. [34]. We just depict 
the 5 < x'^ < 6 interval because the difference between the interacting and noninteracting case 
augments with redshift. 

IV. INTERACTION TERMS OF THE TYPE Q = Q{pc, Pd) 

In case (ii) -i.e., for Q given by Eq. 1^- analytical solution can be found, at most, for one 
of the conservation equations, either (jl]) or ([5]), and this only if = 0, or = 0. Hence, 
at least one these equations must be solved numerically whereby no analytical expression 
exists for the Hubble factor. This limits the usefulness of H{z) to set diagnostics on the 
nature of DE. 

Let us consider first the case F^ = and F^^ > 0. In this instance, the currently phase of 
DE domination would be a transient one, the smaller F^., the longer the DE domination. 

Left and right panels of panel of Fig. O depict the evolution of E"^ and its second derivative 
in terms of for F^ = 0.1 i^o- 

As it can be seen, (PE'^{x)/d[x^)'^ presents a very fast increase (decrease) for phantom and 
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d'E'/d(x')' 
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0.00 




FIG. 5: Dotted line: plot of the second derivative of E'^ for the best fit values, w 
e = 0.0006, with 0^0 = 0.230, and ^ao = 0.725, of model IV, Q = 3He{pc + pd), of Ref. 
line: same as before but with e ~ (no interaction). 



-1.03, 



34(. Solid 



DCC (for quintessence) at small redshifts (z ^ 1). This feature is absent in noninteracting 
models as well as in models whose interaction term is of the general type considered in the 
previous section. 

As a second example, let us consider the interaction term given by Eq. ([6]) with Tc = = 
0.1 Hq. As can be seen in Fig. [71 the behavior of the second derivative, cP (x) / d{x^Y , at 
small redshifts is very similar to the previous case whence we cannot tell one interaction from 
the other. Accordingly, an initial steep slope followed by an abrupt change to a much milder 
slope appears to be a general characteristic of interacting models in which Q is proportional 
to the densities of DM or DE or to a combination of them. We hope this will serve to 
discriminate this set of interactions from the rest. 
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FIG. 6: Left panel: Graphs of E'^ vs. assuming Q = ST^Pfi with = 0.1 Hq for three values 
of w. Right panel: Graphs of the second derivative of E'^ with respect to for the same three 
values of w. In drawing the graphs we took Jl^g = 0.04, 0,cO = 0.26, and ^l^o = 0.70. 




FIG. 7: Left panel: Graphs of E'^ vs. assuming Q = 3{TcPc + ^dPd) with Fc = F^ = 0.1 -ffo 
for three values of w. Right panel: Graphs of the second derivative of E"^ with respect to x^ for 
the same three values of w. In drawing the graphs we took Vt^Q = 0.04, $7c0 = 0.26, VtdQ = 0.70. 



V. THE STATEFINDER PARAMETERS 



Years ago Sahni et al. 



38| and Alam et al. 39|-based on the dimensionless, geometrical 



pair of parameters R and S- gave a somewhat similar diagnostic of dark energy. These are 
defined in terms of the third temporal derivative of the scale factor, namely, R =d / (aH^) 
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and S = (-R — l)/(3 (g — |)), where q stands for the deceleration parameter. The usefulness 
of the said pair stems from the fact that different models are characterized by different 
graphs in the (5, R) plane; in particular, the ACDM model is associated to the single point 
(0, 1). So, they may help elucidate the right DE candidate. 

When these parameters are calculated for cosmological models in which DE and DM 
interact with each other it is seen that the interaction first appears in the third derivative 

n 

of the scale factor [4(| . In particular, for Q = 3e H pc one has 

-R = 1 H W(l + wild) t; 7 TT^ • (31) 

More specifically if one wished to resort to the present value of R to discriminate the 
ACDM model from the DCC model {w = —1, Q = 3eH pc) of sub-section III. A, one would 
find -Ro = 1 — I e ^co- 

If, instead, the interaction is governed by Eq. (j6]) one obtains 

-R = 1 H W(l + wild) t; t; — rr^ . (32) 

2 ^ ^ vib + + {i + w)nd ^ ^ 

This suggests that once the Q parameters come to be accurately determined at different 
redshifts, the statefinder parameters will not only tell whether the interaction really exists 
but also will help discern the expression for Q. However, for the time being, the wide 
uncertainties about these seriously restrict the utility of the R, S pair in this respect. 

Altogether, the diagnostics founded on the and dE'^/dx^ functions appear more ad- 
vantageous since the interaction parameters (i.e., the e and Fj quantities), enter the Hubble 
factor upon integration (analytical or numerical) of the conservation equations ([3])- ([5]). 



VI. DISCUSSION 



Interacting models of dark energy are well motivated, substantially alleviate the coinci- 
dence problem, and show compatibility with observation. Thus the question arises, "can 
they be discriminated from noninteracting models?" For the time being, the answer is not 
given our still imperfect knowledge of basic cosmological parameters, such as the Hubble 
factor. The latter is key in constraining models with observational data from baryon acous- 
tic oscillations, the R-shift parameter of cosmic background radiation, the matter growth 
parameter, etc. Nevertheless, the situation is to improve greatly in the not far future thanks 
to a variety of upcoming and planed experiments. 
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To extract the invaluable information about the nature of dark energy encrypted in the 
history of the Hubble factor is does not suffice to have abundant, accurate data of H{z). One 
must first set up simple and practical criteria to be used once these data become available. 
We proposed several criteria, of notable mathematical simplicity, based on an accurate 
knowledge of the said history -something we reasonably hope to have at our disposal not 
very soon but, at least, in the foreseeable future. 

Specifically if (P E"^ / d{x^Y results positive (negative), then the dark energy is of phantom 
type (either quintessence or a decaying cosmological constant). Only if the DE is not- 
decaying cosmological constant, the graph of E'^ vs. will yield a straight line. 

To discern whether DE is interacting (with Q = Q{Hpc, Hp^)) or not the behavior of 
dE'^{x)/dx^ at high redshifts must be studied. If this derivative tends to the present value of 
the fractional baryon density, Qbo, then there will be grounds to believe that is it interacting 
with the interaction term Q given either by Eq. ([7]) or ( 128|) . If it tends to ^Imo, then we 
may strongly suspect that it is not interacting; and, in particular, if it coincides with Qmo 
independently of redshift, then it will most likely be a conserved cosmological constant -see 
Eq. ([T]). Finally, if the redshift function dE^{x)/dx^ does not tend to any of these two 
quantities, we may conclude that either DE is interacting with Q given by Eq. fl2T]) or 
some other law not considered here, or that in reality w is not a constant but a function 
of redshift. Obviously, this would complicate matters very much because it would call for 
the introduction of further parameters in the analysis -something that we defer to a future 
work. 

If the plot of d"^ E"^ {x) / d{x'^Y ^s. x'^ exhibits a very steep slope at redshifts much smaller 
than unity, we may conclude that the interaction obeys Eq. ([6]) (with Fc and F^ non- 
negative). Its absence excludes such interaction, but it is very hard to specify whether the 
slope is steep enough. On the other hand, since this criterion deals with the slope of the 
second derivative extremely accurate data of H{z) at these redshifts must be employed. 

The statefinder parameters 38|-|40| will be very useful to rule out cosmological models, 
and tell apart different interacting models, but first the histories of VLc{z) and ^d{z) must be 
accurately established in a mo del- independent manner. This looks more challenging than 
obtaining an accurate history of the Hubble factor. 

Since the interaction alters the abundance of dark matter in the past two additional ways 
to ascertain if the DE is interacting or not can be implemented, namely, to examine the 
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behavior of the growth factor, and to study the weak lensing power spectrum. If -as we 
assumed throughout- the DE decays into dark matter, then it must have been less dark 
matter in the past than in the corresponding noninteracting model with the same w and, 
consequently, the growth factor must be smaller. On the other hand, this would imply a 
decline of the lensing spectrum (the opposite would hold true if DM decayed into DE). Both 



effects have been considered by the authors of Ref. 4l|. However, they do not take w 



constant but assume the Chevalier-Polarski-Linder ansatz j42| with very specific values for 
the parameters entering that expression. This is why comparison with our work does not 
seem straightforward. At any rate, both studies may be seen as complementary. 

While the criteria introduced in this paper should be helpful in deciphering the nature 
of DE they should be used alongside the analysis and interpretation of data coming from 
a variety of sources: supernovae type la, cosmic microwave background radiation (CMBR), 
matter power spectrum, baryon acoustic oscillations, weak lensing, and evolution of galaxy 
clusters. The latter lies at the core of recent studies about the dynamical equilibrium of the 



clusters which seem to favor a decay of DE into DM over other possible scenarios 43 |. 

In particular, a promising way to test the existence of the interaction is to consider its 
impact on the integrated Sachs- Wolfe effect. Since the former alters the evolution rate of 
gravitational potentials it must affect the frequency shift experienced by CMBR photons 
crossing collapsing structures. Again, studies of this kind hint at the existence of the inter- 



action with DE decaying into DM 4^, |45| 



We restricted ourselves to constant w and vanishing spatial curvature. While this auto- 
matically limits the scope of our work we feel it does not do it seriously, at least no more than 
otherwise since it would introduce additional unknown parameters that would compromise 
the said scope in other ways. 
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